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Abstract 

Two new analytical solutions of self-induction equation, in Rie- 
mannian manifolds are presented. The first represents a twisted mag- 
netic flux tube or flux rope in plasma astrophysics, which shows that 
the depending on rotation of the flow the poloidal field is amplified 
from toroidal field which represents a dynamo. The value of the am- 
plification depends on the Prenet torsion of the magnetic axis of the 
tube. Actually this result illustrates the Zeldovich stretch, twist and 
fold (STF) method to generate dynamos from straight and untwisted 
ropes. Motivated by the fact that this problem was treated using 
a Riemannian geometry of twisted magnetic flux ropes recently de- 
veloped (Phys Plasmas (2006)), we investigated a second dynamo 
solution which is conformally related to the Arnold kinematic fast 
dynamo. In this solution it is shown that the conformal effect on 
the fast dynamo metric only enhances the Zeldovich stretch, and 
therefore a new dynamo solution is obtained. When a conformal 
mapping is performed in Arnold fast dynamo line element a uniform 
stretch is obtained in the original line element. PACS numbers: 

02.40.Hw-Riemannian geometries 
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I Introduction 



Geometrical tools have been used with success [1] in Einstein general rela- 
tivity (GR) have been also used in other important areas of physics, such as 
plasma structures in tokamaks as been clear in the Mikhailovskii [2] book 
to investigate the tearing and other sort of instabilities in confined plasmas 
[2], where the Riemann metric tensor plays a dynamical role interacting with 
the magnetic field through the magnetohydrodynamical equations (MHD). 
Recently Garcia de Andrade [3, 4] has also made use of Riemann metric to 
investigate magnetic flux tubes in superconducting plasmas. Thiffault and 
Boozer [5] following the same reasoning applied the methods of Riemann ge- 
ometry in the context of chaotic flows and fast dynamos. Yet more recently 
Thiffeault [6] investigated the stretching and Riemannian curvature of mate- 
rial lines in chaotic flows as possible dynamos models. An interesting tutorial 
review of chaotic flows and kinematical dynamos has been presented earlier 
by Ott [7]. Also Boozer [8] has obtained a geomagnetic dynamo from conser- 
vation of magnetic helicity. Actually as pointed out by Baily and Childress 
[9] have called the attention to the fact that the focus nowadays the focus of 
kinematic dynamo theory is on the fast dynamos, specially in construction on 
curved Riemannian manifolds to stretch, twist and fold magnetic flows fila- 
ments or tubes to generated dynamo solutions. This method was invented by 
Zeldovich [10]. In this paper taking the advantage of the success of conformal 
Riemannian geometry techniques used to find out new solutions of Einstein 
general relativistic field equations[l], we use this same conformal geometrical 
technique to find new solutions of the incompressible flows in Arnold metric 
[11, 12, 13]. To resume , finding and being able to recognize the existence 
of dynamos or non-dynamos is not only important from physical and mathe- 
matical point of view, but also finding mathematical techniques which allows 
us to obtain new dynamo solution from dynamos or either nondynamo met- 
rics such as the one we consider here as a flavour of the difficult to find a 
dynamo solution in Riemannian manifolds. This also motivates us to show 
that conformally Arnold fast dynamo metric, is also a dynamo solution. Re- 
cently Hanasz and Lesch [14] have used also a conformal Riemannian metric 
in S 3 to investigate the galactic dynamo also using magnetic flux tubes. Also 
recently, Kambe [15] and Hattori and Zeitlin [16], have investigated the rate 
of stretching of Riemannian line elements of imcompressible fluids, in the 
framework of differential geometry of diffeomorphisms, in the spirit of Zel- 
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dovich stretching. Actually they considered the exponential stretching of 
line elements in time, or dynamo action, in the context of negative curva- 
ture in turbulent flows. They also considered the concentration of vortex 
and magnetic flux tube. This provides a strong physical motivation to the 
present investigation. The paper is organized as follows: In section II the 
the dynamo Riemann metric representing flux rope analytical solution of the 
self-induction equation in the case or zero resistivity is presented. In section 
III the dynamo solution in Riemannian conformal geometry is given. In sec- 
tion IV Riemannian curvature of a particular conformal dynamo is computed 
and in section V conclusions are presented. 



II Thin flux rope dynamos in Riemannian 
manifold 

In this section we shall consider generalization of the Riemann metric of 
a stationary twisted magnetic flux tube as considered by Ricca [17] in the 
Riemann manifold to address the nonstationary case where the toroidal and 
poloidal magnetic fields, in principle, may depend on time. With this metric 
at hand , we are to solve analytical the self-induction magnetic flow equation 
to check for the dynamo existence. Let us now start by considering the MHD 
field equations 

V.5 = (1) 
jjjB- Vx[mxB] -rjV 2 B = (2) 

where u is a solenoidal field while i] is the diffusion coefficient. Equation 
(2) represents the self-induction equation. The vectors t and n along with 
binormal vector b together form the Frenet frame which obeys the Frenet- 
Serret equations 

t = KU (3) 

ft = -Kt + rb (4) 

b' = -Tfi (5) 

the dash represents the ordinary derivation with respect to coordinate s, and 
K,(s,t) is the curvature of the curve where k = Br 1 . Here r represents the 
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Frenet torsion. The gradient operator becomes 



V-f— e- — e— (6) 

ds 60 r 80 6r dr 

Now we shall consider the analytical solution of the self-induction magnetic 
equation investigated which represents a non-dynamo thin magnetic flux 
rope. Before the derivation of this result , we would like to point it out 
that it is not trivial, since the Zeldovich antidynamo theorem states that the 
two dimensional magnetic fields do not support dynamo action, and here as 
we shall see bellow, the flux tube axis possesses Frenet curvature as well as 
torsion and this last one cannot take place in planar curves. Let us now 
consider here the metric of magnetic flux tube 

ds 2 = dr 2 + r 2 dO R 2 + K 2 (s)ds 2 (7) 

This is a Riemannian line element 

ds 2 = g ij dx i dx j (8) 

if the tube coordinates are (r, Or, s) [15] where 0{s) = Or — J rds where r is 
the Frenet torsion of the tube axis and K(s) is given by 

K 2 (s) = [1 - rK(s)cosO(s)} 2 (9) 

Since we are considered thin magnetic flux tubes, this expression is K ps 1 
in future computations. Computing the Riemannian Laplacian operator V 2 
in curvilinear coordinates [16] one obtains 

V 2 = ^[y/ggVdjl (10) 

where dj := and g := detgij where gij is the covariant component of the 
Riemann metric of the flux rope. Here, to better compare the dynamo action 
generation of toroidal field from poloidal fields we shall consider the that the 
toroidal component of magnetic field B s (s) is given in the Frenet frame as 

B s = bo(s)t (11) 
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Note also that we have considered that the flux rope magnetic field does not 
depend on the r and 9r coordinates. While the poloidal magnetic field the 
magnetic field here can be expressed as 

B e (t,0) = e pt b 1 e e (12) 

Now let us substitute the definition of the poloidal plus toroidal magnetic 
fields into the self-induction equation, which with the help of the expressions 

ee = —nsinO + bcosO (13) 

dee*e = — n[(l + T~ l K)sin9 + cos6] — b[cos9 + sinO] (14) 

and 

d t ee = oJtQ — d t nsin9 + d t bcos9 (15) 

Considering the equations for the time derivative of the Frenet frame given 
by the hydrodynamical absolute derivative 

X = d t X +[v.V]X (16) 

where X = (t,n,b) represents the Frenet frame into the expressions for the 
total derivative of each Frenet frame vectors 

t = d t t+[Kb- KTU] (17) 

k = KTt (18) 

b = -K't (19) 

one obtains the values of respective partial derivatives of the Frenet frame as 

d t t = -tk[1 - kt- 2 ^]u (20) 

d t n = tk[1 - KT- 2 ?y}t+ ^b (21) 

db = KT- 1 —n (22) 
r 

where we have used the hypothesis that b = or «'(t, s) = 0, which means 
that the curve curvature only depends on time. A simple example from solar 



5 



physics, would be a flux tube curved and with torsion oscillating with fixed 
sun spots. Substitution of these vectorial expressions into expression (15) 
yields 



r 

along with the equation 



d t eg = [ujcosO -16 — bsinfi — r 1 k— ]n+ [kt(1 + r 2 —)sin6]t (23) 

r r r 



^ = Bqttk (24) 

and the fact that = 0, together with the self-induction equation we obtain 
the following system of equations, for a highly conductive fluid as our own 
universe, with resistivity rj — 

d t B + Tv e sin6B = (25) 

stneB e + ^-[l + r- 2 ^] = (26) 
sint) r 

To obtain these last two expressions we assume that vg >>> v s and that 
d s vs = and that the continuity equation 

V.v = (27) 

where we have considered that the flow is imcompressible which is a reason- 
able approximation in plasma physics. This expresion yields 

^ + Ve rTK = (28) 

OS 

Equation (25) can be rewritten as 

[p + rvesinO] = (29) 
which upon substitution on the equation (26) yields 

^sin9[l + r- 2 ^-} = (30) 

which with the assumption that the flux tube has a small twist and vg 2 << 1. 
Now these equations 
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since B s m constantbo by hypothesis, we have that relation (31) tells us 
that the relation z^f > implies that the poloidal field is amplified from 
the constant modulus toroidal field, which is the dynamo rope condition. 
This physical situation happens often in the sun. Thus relation between the 
angular flow speed and the torsion and radius distance shows that there is 

2 

a lower bound for the thin rope dynamo which is given by r > — . The 
divergence-free equations for the magnetic and flow fields, allows us to write 
down the solutions for Bg 

B e = B°e pt -J i cosdde (32) 

which if we recall the definition of the deviation of flat Riemann metric K of 
the tube above , we may express the integral in terms of K(s) as 

B e = B e pt -I^- K ^ de (33) 

which shows that in the very thin flux rope dynamo in this solution the effect 
curvature of the Riemannian tube is minor. Here we have used the fact that 
the external curvature of the rope is given by = - . The solution for 
torsion and velocity flow are essentially analogous. 



Ill Conformal dynamos on manifolds 

Conformal mapping on a Riemannian line element can in general be repre- 
sented by 

ds 2 = e 2X ^\ds ) 2 (34) 

Manifolds related in this manner, are said conformally related. Note also 
that this is intrinsically connected to the stretch part of STF mechanism to 
generate rope dynamos. Actually in the previous section we considered a 
Riemannian metric for twisted flux tube which was stretched solely on the 
ds element along the magnetic axis of the dynamo rope, through the factor 
K(s), of course when the flux rope is thin, the this stretch effect almost vanish 
though twist and fold may still be kept. Therefore strictly speaking, this is 
not a conformal mapping but only a stretch, therefore strictly speaking not all 
stretches are represented by conformal mapping but every conformal metric 
represents stretching in the Riemannian manifold. One of disadvantages of 
conformal stretching is that the stretching is uniform as in the case of Arnold 
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fast dynamo metric. Conformal metric techniques have also been widely used 
as a powerful tool obtain new solutions of the Einstein's field equations of GR 
from known solutions. By analogy, here we are using this method to yield new 
solutions of MHD dynamo from the well-known fast dynamo Arnold solution. 
We shall demonstrate that distinct physical features from the Arnold solution 
maybe obtained. Before that we just very briefly review the Arnold solution. 
The Arnold metric line element can be defined as [11] 

ds 2 = e~ 2Xz dp 2 + e 2Xz dq 2 + dz 2 (35) 

which describes a dissipative dynamo model on a 3D Riemannian manifold. 
By dissipative here, we mean that contrary to the previous section, the re- 
sistivity 77 is small but finite. The flow build on a toric space in Cartesian 
coordinates (p, q, z) given by T 2 x [0, 1] of the two dimensional torus. The co- 
ordinates p and q are build as the eigenvector directions of the toric cat map 
in 1Z 3 which possesses eigenvalues as Xi — ( - 3+ 2 v ^- > > 1 and X2 = < 1 

respectively. Note for example that if perform a simple constant conformal 
mapping such that 

ds 2 = e 2(Xz+x ^dp 2 + e'^'^dq 2 + e 2X «dz 2 (36) 

which represents a simple global translation and is not changed at every 
point in the manifold. Let us now recall the Arnold et al [13] definition of a 
orthogonal basis in the Riemannian manifold M 3 

e p = e A2 |- (37) 



dp 

,d_ 
dq 



e q = e- A2 |- (38) 



(39) 



d 
dq 

Assume a magnetic vector field BonM 

B = B p e p + B q e q + B z e z (40) 
The vector analysis formulas in this frame are 

Vf=[e Xz d p f,e- Xz d q f,d z f] (41) 
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where f is the map function / : 7Z 3 — > 71. The Laplacian is given by 

A/ = V 2 / = [e 2Xz d 2 f + e' 2Xz d 2 f + d z 2 f] (42) 
while the divergence is given by 

V.B = divB = div[B p e p + B q e q + B z e z ] = [e Xz d p B p + e~ Xz d q B q + d z B z ] (43) 
In particular one may write 

dive p = dive q = dive z = (44) 
in turn the curl is written as 

curlB = curl[B p e p + B q e q + B z e z ] (45) 

where 

curl p B = e~ Xz (d q B z - d z (e Xz B q )) (46) 

curl q B = -e Xz (d p B z - d z ( e - Xz B p )) (47) 

curl z B = e Xz d p B q - e' Xz d q B p (48) 

and 

curle p = —\e q (49) 

curle q = —Xe p (50) 

curle z = (51) 
The Laplacian operators of the frame basis are 

Ae p = —curlcurle p = —X 2 e p (52) 

Ae q = —curlcurle q = —X 2 e q (53) 

Ae z = (54) 

from these expressions Arnold et al [13] were able to build the self-induced 
equation in this Riemannian manifold as 

d t B p + vd z B p = -XvB p + rj[A - X 2 ]B P - 2Xe Xz d p B z (55) 

d t B q + vd z B q = +XvB q + V [A - X 2 ]B P - 2Xe~ Xz d q B z (56) 



d t B z + vd z B z = r)[A - 2Xd z ]B z (57) 

Decomposing the magnetic field on a Fourier series, Arnold et al were able 
to yield the following solution 

b(p, q, z.t) = e Xvt b(p, q,z- vt, 0) (58) 

where B(x,y, z,t) = b(p,q,z,t) and the fast dynamo limit rj = was used. 
Now with these formulas , we are able to compute the solution of the self- 
induced magnetic equation in the background of conformal Riemannian line 
element 

ds 2 = Sl(z) [e- 2Xz dp 2 + e 2Xz dq 2 + dz 2 ] (59) 

The reason for using the general conformal stretching factor Q(z) instead of 
the previous exponential stretching , is to show that the dynamo obtained is 
not only due to the exponential stretching but conformal dynamos, allow for 
the existence of more general conformal stretching. A far obvious, though 
important observation here is the fact that from the equations bellow, we 
recover the Arnold et al [13] if we simply make the conformal factor Q := 1. 
Denoting the dual one form for Arnold basis as 



, = e~ Xz dp 


(60) 


> q = e Xz dq 


(61) 


(fi z = dz 


(62) 



one obtains the Arnold fast dynamo metric as 

ds 2 = <f) 2 + (j) 2 + <p 2 (63) 
Thus the conformal one form dual basis can be expressed as 

<P P C = &<P P (64) 

4> q ° = &4> q (65) 
2 C = &<p z (66) 
On the other hand the vector field basis in conformal metric becomes 

e p = V<-^e Xz — (67) 
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e q = ^e~ Xz §- q (68) 

e 2 = fi-^ (69) 
Let us now repeat some of the fundamental vector analysis relations above 

— * 

in the conformal geometry. The first is the Laplacian of B 

A C B = ft" 1 AS - -tt- 2 [d z tt]d z B (70) 

A fundamental change here in the conformal stretching in Riemannian ge- 
ometrical dynamos, is that of the velocity flow. In Arnold fast dynamo 
example, the flow is a very simple one which is given by (0,0, v) where v is 
constant. Here the dynamo flow is effectively in the dynamo equation by a 
term of the form 

{v.V)B p = n^vdzBp (71) 

which shows that a nonconstant effective velocity flow such as v e ff = Q~ 1 v 
would act in conformal dynamos with respect to the previous Arnold exam- 
ple. The other fundamental component of the curl[vxB] given by 

(B.V)v = n- 2 v[d z n}e z (72) 

With these expressions from conformal geometry in hand, we are now able 
to express the Arnold et al dynamo equations are 

d t B + (v.Vc)B = (B.V c )v + r/A c S (73) 

In terms of components this conformal self-induced equation in the Rieman- 
nian manifold can be expressed as 

dtBp+Q^vdzBp = -Xn^vBp+TiKA-X^Bp^Xe^dpB^-^n-^dztyidzBp+XBp) 

(74) 

The equation for q component can be obtained from p one by simply per- 
forming the substitution A — > —A. The expression for component-z is 

aA + fi -^A = ,[A-2AS s -i^8J B , (75) 
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Decomposing again the magnetic field on a Fourier series now in conformal 
geometry, yields the following solution 

b(p, q, z.t) = e xn ~ lvt b(p, q,z- vt, 0) = e A ^*6(p, q,z- vt, 0) (76) 

where is the conformal Riemannian fast dynamo solution. As given explicitly 
in this solution the basic effect of the conformal geometry in fast dynamos is 
on the speed of dynamo which is an important physical effect. 



IV Riemann curvature of conformal dynamos 

The important role of negative curvature of geodesic flows in dynamos have 
been investigated by Anosov [19]. These Anosov flows, though somewhat ar- 
tificial, provide excelent examples for numerical computation experiments in 
fast dynamos [20]. Within this motivation we include here a simple example 
of conformal spatially stretching, where Q := e 1Xz . This conformal stretching 
applied to Arnold metric yields the conformal metric as 

ds 2 = dp 2 + e 4Xz dq 2 + e Xz dz 2 (77) 

or in terms of the frame basis form of (i = 1, 2, 3) is 

ds 2 = (u p ) 2 + (u«) 2 + (u z ) 2 (78) 

The basis form are write as 

u p = dp (79) 

u q = e Xz dq (80) 



and 



A . 



co z = e i z dq (81) 
By applying the exterior differentiation in this basis form one obtains 

du p = (82) 

du z = (83) 



and 



A 



duji = \e~* z uj z ALU q (84) 
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Substitution of these expressions into the first Cartan structure equations 
one obtains 

TP = = u p q Auj q + u p z Au z (85) 



T q — — \e-2 z u z Au q + u q p Auj p + u q z Ato z (86) 
and 

T z = = uj z p Auj p + u z q Au q (87) 

where T* are the Cartan torsion 2-form which vanishes identically on a Rie- 
mannian manifold. From these expressions one is able to compute the con- 
nection forms which yields 

co p q = -auj p (88) 



. A 



and 



oj q z = \ e -* z uj q (89) 



u)% = (3u p (90) 



where a and j3 are constants. Substitution of these connection form into the 
second Cartan equation 

IV j = R^m^Au 1 = dJj + u)\NJj (91) 

where R % j is the Riemann curvature 2-form. After some algebra we obtain the 
following components of Riemann curvature for the conformal antidynamo 

R% pq = Xe~^ z (92) 

R\ qz = \\ 2 e~ Xz (93) 

and finally 

RP zm = -a\e-$ z (94) 

We note that only component to which we can say is positive is RP zqz which 
turns the flow stable in this q-z surface. This component also dissipates away 
when z increases without bounds, the same happens with the other curvature 
components [21]. 
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V Conclusions 



In conclusion, we have used a well-known technique to find solutions of Ein- 
stein's field equations of gravity namely the conformal related spacetime 
metrics to find a new anti-dynamo solution in MHD three-dimensional Rie- 
mannian nonplanar flows. Examination of the Riemann curvature [21] com- 
ponents enable one to analyse the stretch and compression of the dynamo 
flow. New conformal fast dynamo metric are obtained from the conformally 
self-induced equation.lt is shown that in the effect of conformal mapping 
in Riemannian dynamo flow is to change the fast dynamo speed. Future 
perspectives includes the investigation of homological obstructions in the 
conformal geodesic flows from Anosov flows generalizing the investigation of 
Vishik [22] and Friedlander and Vishik [23]. 
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